Photon generation from vacuum in non-stationary circuit QED 
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We study theoretically the non-stationary circuit QED system in which the artificial atom transi- 
tion frequency has a small periodic modulation in time, prescribed externally. We show that, in the 
dispersive regime, when the modulation periodicity meets the 'resonances', the dynamics may be 
described by the dynamical Casimir effect, Jaynes-Cummings or Anti-Jaynes-Cummings effective 
Hamiltonians. In the resonant atom-cavity regime, under the modulation 'resonance' the dynamics 
resembles the behavior of the dynamical Casimir effect in a vibrating cavity containing a resonant 
two-level atom, and entangled states with two photons can be created from vacuum. Thus, an 
analog of the dynamical Casimir effect may be simulated in circuit QED, and several photons, as 
well as entangled states, can be generated from vacuum due to the anti-rotating term in the Rabi 
Hamiltonian. 

PACS numbers: 42.50.Pq, 37.30.-H, 32.80.Qk 
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INTRODUCTION 



Over the last few decades nonstationary processes in 
cavities have received considerable attention. One such 
process is the nonstationary or dynamical Casimir Effect 
(DCE) - in particular, the creation of photons from vac- 
uum, or another initial field state, in cavity whosegeom- 
etry 0, II H H @ or material properties @, i, i, Hi, 
have a periodic time dependence, with modulation fre- 
quency equal to twice the unperturbed field eigenfrc- 
quency. Nowadays, DCE in cavities is a well studied 
problem, with a variety of theoretical predictions con- 
cerning the number and the statistics of created photons, 
as well as the influence of detuning, dissipation, boundary 
conditions, g eometr y and non-periodicity of the modula- 
tion (see [lCj, HH, G3 for an extensive list of references). 
To date, DCE has not been observed in laboratory, how- 
ever several concrete proposals have appeared over the 

last years Q H IS [H GJ, 01 , with some of them 
being currently implemented experimentally [20]. 

The interest in nonstationary processes in cavities 
reappeared over the last 5 years due to the progress in the 
field of Cavity Quantum Electrodynamics (cavity QED 
|2l| ) in the condensed matter systems, e.g. semiconduc- 
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tor quantum dots [22j, , polar molecules 
perconducting circuits [IE [2(| [27| coupled to resonators, 
the latest architecture known as circuit QED [26|, [28| . In 
cavity QED, the effective two- level atom is coupled to 
the field inside the resonator via the dipole interaction, 
allowing for observation of the light-matter interaction 
at the level of single photons and single atoms. A new 
ingredient in the solid state cavity QED is the possibil- 
ity of engineering and manipulating the properties of the 
artificial atom and the resonator [26|, |29|, |3fJ , as well as 
the interaction strength between them, either during the 
fabrication or in situ. 



Recently, the strong resonant and the strong dispersive 
coupling limits between the artificial atom and a single 
cavity mode were observed experimentally in circuit QED 
[3l| and other solid state architectures [H [23j]. More- 
over, the single photon source [32| . s ing le artificial- atom 
maser (33[, multiphoton Fock states [34| and interaction 
between two artificial atoms (aubits) [30l. |35| were im- 
plemented experimentally, among many other important 
achievements. Besides, the circuit QED architecture ben- 
efits from robust read-out schemes of the atomic internal 
state and the resonator field state [H, HE [H [H, 
relatively low dissipative losses [3lj . state preparation 
techniques [40| and real time manipulation of the atomic 
transition frequency via electric and magnetic fields 
[IE HE HI1 or non-resonant microwave fields [IE [39j] . 

Harnessing the tunability of the atomic transition fre- 
quency, the realization of the Landau-Zener sweeps, with 
the atomic frequency increasing linearly in time, was pro- 
posed in circuit QED 4l| , allowing for generation of sin- 
gle photons and entangled states. In a similar direction, 
in [ll EE S3 the implementation of DCE in semicon- 
ductor cavity QED was considered, using the periodic 
time-dependence of the atom-cavity coupling parameter 
(the vacuum Rabi frequency) . It was shown that there is 
a substantial photon production from vacuum when the 
vacuum Rabi frequency is modulated in time with cer- 
tain 'resonant' frequencies [l7| . A preliminary theoretical 
study of the feasibility of realizing the DCE with a quan- 
tum flux qubit in superconducting quantum nanocircuits, 
as well as the detection of the generated photons, was re- 
ported in . On the other hand, the possibility of con- 
trolling the atomic frequency and detecting the atomic 
internal state is currently being used to couple/decouple 
one or several qubits to/from the cavity mode in order 
to implement quantum logic operations [H HE H^] . 

Here we study the nonstationary cavity QED archi- 
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tecture, in which a single cavity mode is coupled to a 
single artificial atom whose transition frequency has a 
small periodic modulation in time prescribed externally. 
Such a control over the transition frequency, with com- 
patible modulation periodicity, may be achieved in circuit 
QED with present or near- future technology [13, HH [3!| . 
We show that, in the dispersive regime, under the 'res- 
onance' conditions one obtains completely different ef- 
fective regimes for the system dynamics, which may be 
approximately described by the Anti-Jaynes-Cummings 
(AJC), Jaynes-Cummings (JC) or the dynamical Casimir 
effect (DCE) Hamiltonians with adjustable parameters. 
Moreover, in the resonant atom-cavity regime, the dy- 
namics resembles the behavior of the DCE in a cavity 
with oscillating boundaries containing a resonant two- 
level atom (detector) [44[, and entangled states with up 
to two photons can be generated from vacuum. 

Thus, we demonstrate the possibility of simulating the 
DCE in circuit QED using a single non-stationary atom, 
instead of a macroscopic dielectric medium as in [6[. As 
applications, it may be possible to create excitations, ei- 
ther photonic or atomic, from the initial vacuum state 
1 3,0), generate nonclassical states of light and realize 
transitions between the states {\g,m}, |e,m± 1)} in the 
dispersive regime. Here \g) and |e) stand for the atomic 
ground and excited states, respectively, and |m) for the 
Fock state of the cavity field. A related problem was 
recently studied in (4f|, where it was suggested that las- 
ing behavior and the creation of a highly non-thermal 
population of the oscillator, as well as the cooling, could 
be implemented using an analogous scheme in the near 
future. 



II. NONSTATIONARY CIRCUIT QED 

We assume that the atomic transition frequency (t) 
may be described as the sum of two terms. The first term 
fio describes the bare atomic frequency and the second 
term represents a small modulation amplitude e <C f^o 
multiplied by a periodic function of time ft prescribed 
externally 



H (t) 



where 



= at, 



|e)(e| — | g) (g | and we assume 



H = 1. In the stationary case, s = 0, one may perform 
the Rotating Wave Approximation (RWA) and obtain 
the standard Jaynes-Cummings (JC) Hamiltonian [4(| . 
which has been verified in several experiments over the 
last few years [H, HH, H3, EH, HU . However, in the non- 
stationary case, as well as under strong dephasing noise 
[50l ]. the anti-rotating term (acr_ + aV + ) cannot always 
be eliminated. Moreover, it is responsible for producing 
an analog of the DCE and creating photonic and atomic 
excitations from vacuum under modulation 'resonance' 
conditions, as shown below. 

In the interaction picture with respect to Hq (t) the 
interaction Hamiltonian reads 



Hi = 30 (e l 



(3) 



where h.c. stands for the Hermitian conjugate and S± = 
J dr [il (t) ± ui] . All the information about the dynam- 
ics of the system is contained in the time-dependent co- 
efficients exp(zS±), which may be significantly simpli- 
fied by adjusting the modulation frequency r\ in order to 
achieve the 'resonances'. We have explicitly 



3 e iH± = 3e iA±t 



En 



l=Q L ' k=l 



- (^ke- lknt - A* k e ikr,t ) 



(4) 

where we defined a complex coupling constant g = 
go exp \i (e/rf) Y^kLi k _1 Sfc] an d parameters 



At. 



Cfc + isk 
2k : 



A, 



■ eco ± u>. 



III. AJC AND JC RESONANCES 



The 'Anti-Jaynes-Cummings' (AJC) resonance occurs 



for 



n(i) = fi +e/t, ft = y~) (sk sin k-qt + c k cos krjt) . (1) 

Here r\ is the modulation frequency and {sk,Ck\ form a 
set of coefficients describing an arbitrary periodic time- 
dependence of ft- We suppose that the cavity frequency 
u) and the atom-cavity coupling parameter 30 are con- 
stant, so at the 'sweet spot' in solid state cavity QED 
the system is described by the Rabi Hamiltonian (RH) 

M " 



H = Hq (t) + 30 (a + a f ) (ct+ + <j_) , 



(2) 



where a (a^) is the cavity annihilation (creation) opera- 
tor, <t + = \e)(g\ and er_ = |<7)(e|. The free Hamiltonian 



V = Vajc = a +~^, 

where |£| <C t] is a small 'resonance shift'. Assuming a 
reasonable experimental condition ejr\ <^ 1 we expand 
(j4]) to the first order in ejr\ and make the RWA in ((3|), 
obtaining 

Hi ~ 3 (6»e^*aV + + e iA - t aa+) + h.c, (5) 

where the AJC dimensionless coupling is 

= A ie /?7. (6) 

In the resonant regime, |A_|/3o <SC 1, we apply the 
method of slowly varying amplitudes to the Hamilto- 
nian ([5]) , repeating the procedure employed originally for 



3 



studying the photon generation from vacuum due to the 
DCE in a vibrating cavity containing a stationary two- 
level atom [3]. We find that for the initial state \g, 0) the 
photon generation occurs for two values of the resonance 
shift 

C = e±- A_/2±V2. 9o 

and one gets the following non-zero probabilities P x ,m> 
with x = {e,g} denoting the atomic state and m the 
photon number 



P gfi « cos 2 (xt) , Pes w sin 2 (y + q) sin 2 ( X t) , 
' 3) 2 w cos 2 [y - q) sin 2 (xt) . (7) 



P 



Here 



X ~ do \0\ sin (y + q) , tany : 



2\/2go + A- 
2V2 5 o ~ A. 



1/2 



and q = (7r/2) for £_ (£+)■ Thus, in the resonant 
regime, when the atomic transition frequency is modu- 
lated with the frequency A + — £±, a superposition of 
states |<7,0), |e, 1) and |y,2) is created from the initial 
vacuum state \g,Q), and the probability of exciting the 
atom is limited by sin 2 y sa 1/2. 

We illustrate this behavior in Fig. la, where we show 
the exact dynamics of P g fi, Pe,i and P g ,2 vs. time (t) 
for the AJC resonance, ft = smr] A j C t with £ = us- 
ing the parameters go/u> — 4 • 10 -2 , A_ = go/ 10 and 
e/oj = 10 _1 . This dynamics resembles the one occurring 
in the context of DCE [44| . where a resonant (station- 
ary) two-level atom or detector is fixed inside a cavity 
whose boundary is oscillating with the frequency close to 
2u). In both cases not more than two photons can be cre- 
ated from the vacuum state |<7,0) and the probability of 
exciting the atom is limited by the value 1/2. This simi- 
larity is not surprising, since in this case the modulation 
frequency A + w 2ui, and the atom plays the role of the 
two-level detector and the cavity modulating mechanism 
(via the atom-cavity coupling) at the same time. 

In the dispersive regime, *7o \J ( n ) I |A_| <C 1, where (n) 
is the mean photon number, the Hamiltonian ([5]) may be 
approximated by [46] 



H 



(i) 



9e^Vcr + + h.c.) +5{n + 1/2) a z 



where 



S = gl/^ 



is the dispersive shift. In a rotating frame we obtain the 
effective AJC Hamiltonian 

H AJC ~ [£ + 6(1 + 2n)] y + {gOa'a + h.c) (8) 

By adjusting the resonance shift £ [52I ] in order to make 
|£ + <5 (1 + 2n)| small compared to |<70|, one obtains the 
resonant AJC Hamiltonian. From the physical point of 
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FIG. 1: a) Exact dynamics of P g ,o, Pe,i and P 9i 2 vs. time 
in the resonant atom-cavity regime under the AJC resonance, 
b) Exact dynamics of (n) , P g and P e in the dispersive regime 
for the AJC resonance and c) JC resonance. 



view, the external modulation supplies the energy lu + firj 
necessary to create one photon and one atomic excitation 
simultaneously. Thus, one can create the superposition of 
states |e, 1) and \g, 0) starting from the initial state \g, 0), 
as illustrated in Fig. lb, where we show the exact dynam- 
ics of (n) and P g for the AJC resonance, f t = smrj A j C t 
with £ = —2(5, using the experimental circuit QED pa- 
rameters Qq/uj — 1.4, go/u) = 2 • 10~ 2 and assuming 
e/w = 0.2. 

Also in the dispersive regime, the 'Jaynes-Cummings' 
(JC) resonance occurs for 

r, = rij C = |A_| - £. 

For positive A_ we get the effective JC Hamiltonian 

[£ + 5(l + 2n)]^ + (g0 aa++h.c). (9) 



H 



JC 



If A_ is negative, we obtain the same effective Hamil- 
tonian upon replacements 8^—8* and £ — > — £ in Eq. 
©. Thus, by employing the JC resonance and adjust- 
ing the resonance shift £, one may couple the subspaces 
{|<7,m), |e, m — 1)} when the atom and the field are far 
off- resonant, since the external modulation supplies the 
energy difference \co — Oo| necessary to couple the atom 
and the cavity field. This behavior is illustrated in Fig. 
lc, where we show the exact dynamics of (n) and P e for 
the JC resonance, ft = smr]j C t with £ = —26, for the ini- 
tial state \g, 1) and the parameters of Fig. lb. Moreover, 
one could engineer entangled states with several photons 
from the initial vacuum state \g, 0) by alternating be- 
tween the AJC and JC resonances and controlling the 
time interval and the resonance shift of each resonance. 
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IV. DCE RESONANCE 

In the dispersive regime, the 'dynamical Casimir effect' 
(DCE) resonance occurs for 

V = Vdce = 2w - 2£. 

Performing the RWA in the interaction Hamiltonian 
(|3|), in a rotating frame we obtain the time- independent 
Hamiltonian consisting of the JC Hamiltonian plus the 
AJC term multiplied by the adjustable coupling g9 

ffj 1 ' ~ £n+ A ~ 2 + - a z + (gaa + + g6a)(T + + h.c.) . (10) 

We may obtain an effective Hamiltonian by applying a 
sequence of small unitary transformations [46[ on l|10p 
and performing the Hausdorff expansion after each step. 
Assuming that 6 ~ 0(\g/A_\) we apply the 'rotating' 
unitary transformation 

U r = exp [(gaa+ — h.c.) /A_] 

followed by the 'anti-rotating' one 

U a = exp [(g9a^a+ - h.c.) / (A_ + 20] 

to obtain the effective Hamiltonian H e ff = 
U a U r Hj U^U^, which in a rotating frame reads 

H ef f ~ {£ + 6<j z )n + S<T z (9*a 2 + h.c.) (11) 

- |^ {ge lA - t ana + + h.c.) + 0(\g/A^\ 3 ). 

The first two terms of the effective Hamiltonian (fTTj) 
form the DCE part and the remaining terms represent 
the corrections, whose leading term (oscillating with high 
frequency ~ |A_|) describes the nonresonant photon ab- 
sorption by the atom. These corrections become relevant 
when the third term becomes large, so for initial times 
(roughly for go •*/ (n) / |A_| <C 1) their contribution is rel- 
atively small and a z becomes approximately a constant. 
If the atom is initially in the ground state, the Eq. (fTTj) 
becomes the DCE Hamiltonian 

H DC E^({-5)n~6(0*a 2 + h.c). (12) 

For the atom initially in the excited state, a similar effec- 
tive Hamiltonian is obtained under substitution S — ► —8. 
Therefore, by adjusting the frequency shift to £ = ±<5, 
depending on the initial atomic state [53j . we have pho- 
ton pairs creation from vacuum and field amplification 
due to an analog of the DCE. 

Here the DCE is simulated by the atomic transition 
frequency modulation through the atom-cavity coupling. 
However, the photon generation process is not steady be- 
cause after several photons have been created the third 
and further terms in (|lip become important, and the pho- 
ton generation is interrupted. Nevertheless, the Hamil- 
tonian (jlip shows that it is possible to simulate DCE 




FIG. 2: Exact dynamics of (n) and P e vs. time in the dis- 
persive regime for the DCE resonance. {n) DCE is the mean 
photon number for the DCE Hamiltonian (|12p . 

and generate several photons from vacuum using a sin- 
gle artificial atom. This is illustrated in Fig. 2, where 
we show the exact dynamics of (n) and P e vs. time for 
the DCE resonance, ft — sin rj DCE t with £ = 5, for 
the initial state |<?,0) and the parameters £Iq/lu = 1.4, 
go/oj = 2 ■ 10~ 2 , e/uj = 0.4. We also show the curve 
( n ) dce = s i nrj2 (2 \S6\ t), which gives the expected mean 
photon number for the DCE Hamiltonian (|12p . demon- 
strating that for initial times the exact dynamics can 
be described by the DCE Hamiltonian. However, after 
the creation of a few photons the atom acquires a finite 
probability of being excited, and the dynamics starts to 
deviate from the DCE Hamiltonian. The photon gener- 
ation is interrupted and restarts again as time goes on, 
respecting the limit 9oVw/ |A-| <C 1, and the behavior 
of P e resembles the one of (n) . 

This phenomenon may be qualitatively understood as 
follows. In the dispersive regime the atom acts as an 
effective non- linear capacitance [36(, pulling the cavity 
frequency to 

Q (i) w w + . ° z9 } , « (w + a z 5) - a z S^-f t . 

A_ + ejt A_ 

Consequently, one expects that the periodic modulation 
of ft with the modulation frequency close to rj sa 2 (u> ± S) 
would lead to DCE [ll|, [44|, for which the photons are 
generated as long as the modulation is present. The 
energy 2uj necessary to create pairs of photons is pro- 
vided through the atomic frequency modulation and the 
resonance shift £ must be adjusted [!U m order to get 
a constructive interference on the cavity field [TTI . 
However, as time goes on the atom gets entangled with 
the field, acquiring a finite probability of being excited 
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through photon absorption, and the photon generation 
cannot continue steadily due to the loss of constructive 
interference. This is different from the usual DCE sit- 
uation, where the properties of the macroscopic linear, 
lossless and nondispersive dielectric medium inside the 
cavity are modulated 0, [1, @, [|[ , and the field does not 
get entangled with individual atoms. 

V. DISCUSSION AND CONCLUSIONS 

In the previous sections we have considered only the 
first order resonances. In general, the K-th order reso- 
nances occur for an integer K when 



where r\ i stands for the AJC, JC and DCE resonances. 
In this case one recovers the previous result upon substi- 
tutions 5 — ► 5k ~ 5 + gi/A+ ■ ■ ■ and 9 — ► 9 k, as follows 
from the expressions and ([?]). Now the effective dis- 
persive shift 5k contains the contribution of many terms, 
among them the Bloch-Ziegert [46| shift g^/A + and pow- 
ers of {e/r/). 9k contains contributions due to the non- 
harmonic shape of the pulse, A fc>1 , as well as due to the 
powers of {e/r}) 

9 K = Ak e/v + --- + (Ai e/nf /K\ + ---. 

However, to employ the higher order resonances the ef- 
fective dispersive shift 5k should be carefully evaluated, 
otherwise there is a risk of missing the exact modula- 
tion resonance, since \9k\ becomes smaller and there is 
less freedom in committing small errors in the resonance 
shift £. 

Our results may be easily translated to the situation 
where f2 {t) — f2o is constant and go {t) = go + eft has 
a periodic time-dependence. In this case the interaction 
Hamiltonian is 

Hi = go (*) (e^°-^W+ + e l (° 0+ ^ t a t (j + + /i.e.) . 

(13) 

If one expands exp{i^±) in Eq. (H} to the first order 
in e/r}, the Hamiltonian ([3]) becomes equivalent to the 
Hamiltonian (fTB")) . so the results obtained above for {t) 
also hold for go (t) after making appropriate substitu- 
tions. The main difference is that in the O {t) case the 
higher order resonances occur due to the powers of e/r/ 
and non-zero coefficients Ak , while in the go {t) case they 
are due to non-zero coefficients Ak only. Finally, if the 
cavity frequency is modulated periodically, with constant 
fi and go, the AJC and JC resonances also occur, besides 
the well known DCE resonance [iil ]. 

The experimental verification of this scheme seems pos- 
sible in circuit QED architecture with superconducting 
qubits and coplanar waveguide resonators [36| . where 
one may adjust the system parameters in situ via elec- 
tric and magnetic fields, as demonstrated in [261 130. l35j|. 



Moreover, several schemes to read out the ca vity and the 
atomic states are currently available [3l|, [H, [38| and un- 
der investigation [3!| . The main issue would be to mod- 
ulate periodically the atomic transition frequency with a 
stable modulation frequency r/ ~ 10 GHz, what is within 
experimental reach 30]. One could also use this scheme 
to couple M identical qubits (e.g. superconducting 2- 
level atoms [30j, |35[ or a cloud of polar molecules [24[) 
to the same cavity mode and modulate the frequency of 
M atoms simultaneously, since in this case the effective 
coupling increases to \fMg . 

One important point we did not analyze here is the 
dissipation and decoherence of both the artificial atom 
and the cavity due to the noisy solid state environ- 
ment {3^ . l5ll ]. Recent experiments achieved experi- 
mental values {k/uj < 10~ 4 ,7/cj < 10~ 3 ,7 ph /w < 10~ 3 } 
[3ll | . where k is the cavity decay rate, 7 is the atomic de- 
cay rate and 'y ph is the atomic pure dephasing rate. To 
deal with dissipation in a qualitative manner, we com- 
pare the rates of the photon production from vacuum 
for each resonance to the dissipation rates. We take 
the current experimental value of the coupling constant 
go/u ~ 2 • 10~ 2 [31] and assume e/uo ~ A_/w ~ 10 _1 to 
make the estimative. The photon creation rate for the 
first order DCE resonance is roughly \59\ /to ~ 10~ 4 , and 
for the first order AJC resonance \g9\ /u> ~ 10~ 3 . Both 
these values are larger or of the order of magnitude of the 
dissipation/decohcrence rates, so the photon production 
due to modulation of Q {t) seems possible in the future. 

In conclusion, we analyzed the nonstationary circuit 
QED system in which the atomic transition frequency 
has a small periodic modulation in time, prescribed ex- 
ternally. In the dispersive regime, under the modulation 
'resonances' the dynamics can be effectively described 
by the Anti-Jaynes-Cummings, Jaynes-Cummings or the 
dynamical Casimir effect Hamiltonians. Moreover, in the 
resonant atom-cavity regime, under the corresponding 
'resonance', an entangled state with two photons can be 
created from the vacuum state |g,0), analogously to the 
dynamical Casimir effect in a vibrating cavity containing 
a two-level atom. This study illustrates the importance 
of the anti-rotating term in the Rabi Hamiltonian, ne- 
glected in the Jaynes-Cummings model - here this term is 
responsible for photon generation from vacuum and field 
amplification. As applications, this scheme can be used 
to verify photon creation from vacuum in nonstationary 
circuit QED due to an analog of the dynamical Casimir 
effect using a single atom, as well as off-resonant transi- 
tions between the states {\g,m), |e,m± 1)} and genera- 
tion of entangled states with several photons. 
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